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Abstract: We present a novel method for efficient sorting of photons 
prepared in states of orbital angular momentum (OAM) and angular 
position (ANG). A log-polar optical transform is used in combination with 
a holographic beam-splitting method to achieve better mode discrimination 
and reduced cross-talk than reported previously. Simulating this method 
for 7 modes, we have calculated an improved mutual information of 2.43 
bits/photon and 2.29 bits/photon for OAM and ANG modes respectively. 
In addition, we present preliminary results from an experimental imple- 
mentation of this technique. This method is expected to have important 
applications for high-dimensional quantum key distribution systems. 
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Sorting photons according to their transverse spatial mode is an interesting problem and the 
subject of ongoing research Q). Perfect sorting occurs when photons prepared in orthogonal 
input modes are transformed into modes whose intensities do not overlap in space or in time. 
The ability to accurately and efficiently discriminate the transverse modes of individual photons 
is vital for applications seeking to encode quantum information using the spatial degree-of- 
freedom. Sorting methods should aim to minimize detection errors and loss. 

Modes with helical phase fronts and no other azimuthal dependence are of particular inter- 
est because their helical structure is preserved in propagation through cylindrically symmetric 
systems, for example in the free-space propagation between an optical transmitter and receiver 
with circular pupils. It has been long known that modes with azimuthally dependent phase e f 
carry £h of orbital angular momentum (OAM) ||2l[3]. The complex field of such a mode can be 
represented by 

u t (r,<p)=R(r)J t *, (1) 

where R(r) is an arbitrary function of the radial coordinate r and £ is an integer. Owing to their 
high-dimensionality, OAM modes are a prime candidate for free-space quantum key distribu- 
tion (QKD) systems As such, the ability to efficiently sort single photons based on their 
OAM mode number has become the focus of current research. 

Recently, a method for discriminating light beams based on their OAM quantum number I 
has been demonstrated [ 8 9 1 . Two phase-only holograms are used to optically map polar coordi- 
nates (r, (p ) in the input plane to rectilinear coordinates in the output plane (x,y) via the log -polar 
mapping x = a(q> mod 2n) and y = — aln(r/b), where a and b are scaling constants 11 1 011 111 . 
Thus, optical waves with helical phase fronts are transformed into tilted plane waves, which 
can be sorted at the focus of a lens. While this method was substantially more efficient than 
previous methods for sorting OAM modes, it was still limited to an efficiency of approximately 
80% (9). That is to say that for a photon with OAM I, there existed an approximately 20% 
probability of detecting it with OAM m 7^ £, 

In this work, we show that the technique of Berkhout et al. [8] can be combined with a 
holographic beam-splitting technique to sort OAM modes with near unit efficiency. We do this 
by modifying the two phase-only elements used in Ref. [8] appropriately. In addition, we show 
that a similar method can be used for efficient sorting of modes complementary and unbiased 
to the OAM basis. These modes are represented as a superposition of OAM modes 

1 L 

1 ^ - -i2^/(2i+l)_ 



V + 1 ( = _l 



and are referred to as angular modes (ANG) because of the angular localization of their intensity 
patterns. The OAM and ANG bases, by being mutually unbiased with respect to one another, 
guarantee security against eavesdropping in QKD due to the inability of an eavesdropper to 
detect a photon simultaneously in both bases [4, 12 1. Thus, for an OAM-based QKD system, an 
efficient method of sorting photons in the OAM as well as the ANG basis is essential. 

To proceed, we consider in detail the action of the OAM sorter [8 1 on the OAM modes given 
by Eq. (|T). Without loss of generality, we drop the radial dependence and treat the problem 



only in the azimuthal dimension. After applying the log-polar coordinate transformation, these 
fields can be expressed as the truncated plane waves 

(4), (3) 

where rect(x) = 1 for |jc| < 1/2 and otherwise. At the focus of a lens with focal length /, we 
can express the field as 

/x'-AA 

U e (x') = (27ta) sine ( — — J , (4) 

with sinc(x) = sin(7Dc) / (nx) and A = fXj (2na). Thus, an input field with OAM index of I is 
focused to a spot centered about x' = At. Due to the non-zero width of these spots, there is a sig- 
nificant amount of overlap between the neighboring modes in the output intensity pattern. This 
shows that even in principle, the OAM sorter cannot perfectly discriminate between adjacent 
OAM modes (9). 

However, the amount of spatial overlap between the output modes can be reduced. The key 
is to recognize that the blurring of the spots in the detection plane is caused by the use of only 
a single angular cycle of (p (8). To overcome this, we instead periodically map several cycles 
of <p to the points in the output coordinate x. For example, the field at angle <p maps to points 
x = a(q> + 2nm) for all integers m. One method to accomplish this mapping is to coherently 
split the field given in Eq. (0 into an array of N orders spaced by 2na. This results in the field 
U' t (x) = e itx l a re,ct(x/(2naN))/\/N that when focused yields 

~„ , s 2%a . fx!-Al\ 

U> e (x>) = — smc^j. (5) 

Hence, the resulting overlap becomes negligible when N becomes large enough and perfect 
discrimination of OAM fields can be realized. 

Angular (ANG) modes can also be sorted using a similar approach. The coordinate transfor- 
mation acting on mode 9j yields 



©,M = ^^rect(^Wfc^ 
V2L + T \2na) V a 



(6) 



where A e = 2na/{2L+ 1) and 8 L {u) = ^ = _ L e itu is the Dirichlet kernel. The intensity pattern 
has a peak at x = Ag j and its first zero occurs at x = Ag (j ± 1). Thus, as with the OAM modes, 




Fig. 1: (a) ID profile of the fan-out phase element designed to diffract an incoming beam into 9 uni- 
formly spaced copies (b) ID profile of the phase-correcting element designed to correct the relative phase 
introduced by the fan-out element between the 9 copies. 
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Fig. 2: Schematic showing the configuration for enhanced sorting of ANG modes. The OAM sorter shown 
here is from Ref. [8 1. The fan-out phase element is combined with the second element of the OAM sorter 
to create multiple copies of the transformed beam. A phase-correcting element then corrects the relative 
phases introduced between the copies. For enhanced sorting of OAM modes, the fan-out phase element is 
directly combined with the reformatter 1 8] and one phase correcting element is used. 



ANG modes can be sorted using position information. However, the spots corresponding to 
neighboring modes overlap spatially, leading to crosstalk in detection. 

Like the OAM modes, the resolution of the transformed ANG modes is ultimately restricted 
by their limited spectrum. Their spatial frequency content is given by the Fourier transform of 
Eq. © 

& (k) = 2%a y sine (ak - i) e - W(2£+i). (7) 

Analogous to the technique we used with the OAM modes, we can extend the spatial bandwidth 
of the modes in such a way that the performance of the sorting improves. 

Using a lens of focal length /, we focus the transformed angular mode, yielding the field 
©'■ (2nx! /Xf). If we now coherently split this field into an array of N orders spaced by (2L + 
\)Xf/(2nd) and refocus this field using another lens of focal length /, the resulting field will 
be identical to the field in Eq. © but instead the sum will contain N(2L+ 1) terms. When 
N = 2n + 1 is odd, the field is written as 

= , 1 rect (JL) 8 nL ( X -^l) , (8) 
; y/N(2L+l) \2naJ nL \ a J' 

which has a narrower width than the field in Eq. (|6]l by a factor of N. Consequently, arbitrarily 
low crosstalk can be achieved in the sorting of ANG modes by choosing a sufficiently large N. 

The implementation of this technique is straightforward. In addition to optics required for 
the log-polar coordinate mapping, we require an efficient method of splitting the field into N 
copies. The fan-out element introduced in Ref. |fT3l is a phase grating designed to diffract an 
incoming beam into N uniformly spaced orders, each having the same spatial profile and equal 
energy. For perfect beam splitting, an optical element has to transform an incoming plane wave 
into a field distribution given by 

N 2%s m 
U{x,y) = £ A m exp(i<j) m )exp(-i—r—x), (9) 

m 1 

where A m is the amplitude, <p m is the phase, and s m is the angle of propagation of the Af copies. 
The fan-out element is the optimal design in the family of phase-only holograms which can 
approximately achieve this task |fl4l . Generally, the fan-out element will introduce a relative 
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Fig. 3: Simulation results comparing the output from the OAM sorter (broad peaks) with the output 
from the N = 9 fan-out-enhanced OAM sorter (narrow peaks) for 7 input OAM modes. Different colors 
correspond to different input OAM modes. 

phase (j) m between the different copies. This can be removed with a phase-correcting element 
in the Fourier plane of the fan-out element. The multiple copies are then Fourier transformed 
with a lens to a narrower spot than before. Using the specific values of A m and (j> m given in 
Ref. Ifl3ll and 03), we can achieve an efficiency of more than 99% while splitting the beam to 
nine copies. One-dimensional profiles of the fan-out and phase-correcting elements are shown 
in Fig.Q] 

These two phase-only elements are combined with the two elements used for sorting OAM 
modes in Ref. |8|. Then, the first element in our enhanced technique performs a log-polar 
optical transformation, as well as creates multiple copies of the transformed modes. The second 
element is used to correct the phase error introduced by both transformations. As in Ref. JHJ , 
a Fourier transforming lens is used between these two elements. For sorting ANG modes, the 
fan-out element is placed directly in the exit plane of the coordinate mapping optics instead 
(Fig. [2}. The relative phase introduced between the multiple copies is then corrected by a third, 
separate element. For both methods, a final lens is used to refocus the light onto the detector. 
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Fig. 4: Simulation results comparing the output from the ANG sorter (broad peaks) with the output from 
the N = 9 fan-out-enhanced ANG sorter (narrow peaks) for 7 input ANG modes. Different colors corre- 
spond to different ANG modes. 
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Metric\Sorter 


Limit 


OAM 


Fan-Out (OAM) 


ANG 


Fan-Out (ANG) 


Mutual Information 


2.81 


1.75 


2.43 


1.55 


2.29 


Error Percentage 


0% 


20% 


5% 


23% 


7% 



Table 1: Mutual information (in bits/photon) and error percentages calculated for the OAM and ANG 
sorters and the fan-out-enhanced OAM and ANG sorters used for sorting 7 modes. The theoretical limit 
is shown in the second column and is equal to log 2 (7) =2.81. 



Figure [3] shows simulation results comparing the output of the OAM sorter from Ref. [8| 
with the output from the fan-out-enhanced OAM sorter, for different input OAM modes. Figure 
|4] shows similar results for the ANG sorter for different input ANG modes. The change from 
broad, overlapping peaks to narrow peaks with very little overlap clearly shows the dramatic 
improvement in sorting ability. We quantify this improvement in Table 1 by calculating the 
mutual information and error percentage for each method. The error percentage refers to the 
probability of obtaining mode x ^ y in the output if mode y was input. The mutual information 
is calculated by the formula 

I = H(x)-H(x\y) = -^A-log^PO + L^L^v^^vx). (10) 

x x y 

where the summations are performed over the total number of modes in our system. P x is the 
probability that mode x was input, and P yx is the conditional probability that mode y was output 
given that mode x was input. From Table 1 , one can see that the percentage error is substantially 
reduced in both cases and the mutual information is increased from 1.75 to 2.43 bits/photon for 
the OAM sorter and from 1.55 to 2.29 bits/photon for the ANG sorter. While these results show 
a definite improvement from previous sorting methods [ 8 , 9 1 , the fan-out method can reduce the 
error further and bring the mutual information towards the theoretical limit of log 2 (n) . This can 
be accomplished by creating more copies by using two cascaded fan-out elements fl3l . 

We are currently working towards an experimental implementation of this method using two 
phase-only spatial light modulators. Fig. [5] shows preliminary results comparing the output 
from the OAM sorter of Ref. [9] with our enhanced OAM sorter. Here, the fan-out element 
creates N = 3 copies of the incoming beam. As seen above, this technique should allow us to 
realize near-perfect sorting of OAM and ANG modes, which will prove invaluable for quantum 
information protocols employing OAM modes BT fTSI . This work was supported by the DARPA 
InPho program and the Canadian Excellence Research Chair (CERC) program. 
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Fig. 5: Preliminary results from an experimental implementation of our sorting method, (a) Output from 
the OAM sorter |9| for different input OAM modes, (b) Output from the fan-out-enhanced OAM sorter 
proposed in this paper. Here, we used a fan-out element that creates 3 copies of the beam. 





